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We show that, under a reasonable "smoothness" hypothesis, it is possible to extract informations on the 
amplitude and phase of two-body non-leptonic weak decay matrix elements from the study of Euclidean correlation 
functions in lattice QCD. 
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1. INTRODUCTION 

Exclusive non-leptonic weak decays are a fun- 
damental ingredient in the determination of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix el- 
ements and in the study of CP violation in 
the K-, D- and B-meson decays. However, the 
study of these non-leptonic decays involves non- 
perturbative aspects of the strong interactions, a 
challenging task for our still incomplete knowl- 
edge of non-perturbative dynamics. Unfortu- 
nately, a theoretical description of exclusive de- 
cays based on the fundamental theory is not pos- 
sible yet. Over the years, several methods have 
been introduced to estimate the relevant matrix 
elements: vacuum saturation, bag models, quark 
models, QCD sum rules, 1/N C expansion, chiral 
Lagrangians, factorization, etc. However, none of 
these approaches is actually based on first princi- 
ples. 

The lattice approach has been used to obtain 
results based on first principles for a wide set of 
relevant physical quantities such as the hadron 
spectrum, the meson decay constants, the form 
factors entering in semileptonic and radiative de- 
cays, the kaon i3-parameter Bk, etc. However, 
any computation of exclusive non-leptonic de- 
cays performed in the Euclidean space-time suf- 
fers from severe limitations, as shown by Blok and 
Shifman jj| in the context of QCD sum rules, and 
by Maiani and Testa j|] in the framework of lat- 
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ticc QCD. In fact, the activity in this field HQ 
has completely stopped after the publication of 
ref. |. 

In this paper we show that, under quite reason- 
able physical hypotheses, it is possible, at least in 
principle, to extract predictions for the relevant 
matrix elements in numerical simulations of lat- 
tice QCD, in spite of the difficulties due to the 
Maiani-Testa No-Go Theorem (MTNGT). 

The MTNGT states essentially the following: 

• In the calculation of a two- (many-) body 
decay amplitude performed in the Eu- 
clidean space-time, which is the only pos- 
sibility in Monte Carlo simulations, there is 
no distinction between in- and cwf-states. 
As a consequence, the matrix elements that 
one is able to extract are real numbers re- 
sulting from the average of the two cases. 
This jeopardizes the possibility of any re- 
alistic prediction for the matrix elements. 
For example, we know from the measured 
A1/2 and A 3 / 2 amplitudes in D — > Kir de- 
cays that there is a phase difference of about 
80°. 

• Matrix elements are extracted on the lattice 
by studying the time behaviour of appropri- 
ate correlation functions at large time dis- 
tances. Maiani and Testa showed that what 
can be really isolated in this limit are the 
off-shell form factors corresponding to the 
final particles at rest. For kaon decays, we 
can use the chiral theory to extrapolate the 
form factor to the physical point. This is 
certainly not the case for D- and B-meson 



decays. In the latter case it is not possi- 
ble to obtain a realistic prediction for the 
matrix element. 

We will show that both these difficulties can 
be overcome under a "smoothness" hypothesis, 
and that under this hypothesis it is possible, at 
least in principle, to extract the physical matrix 
elements, including the phase due to the strong- 
interaction rescattering of the final states. In the 
interesting case in which Final State Interaction 
(FSI) is dominated by the exchange of a reso- 
nance in the s-channel, it is possible to calculate 
the parameters of the resonance, as explained in 
detail in ref. ||. 

Unfortunately, at present we do not know if the 
method we are proposing can be succesfully ap- 
plied in numerical simulations on currently avail- 
able lattices. 

2. CORRELATION FUNCTIONS IN 
THE EUCLIDEAN SPACE-TIME 

Following ref. Q, we first examine the Eu- 
clidean three-point function G$(t\, t 2 ): 

Gq(ti,t 2 ) = (OlTp^tijn^tiOffCOJllO) 

= (o|n g -(t 1 )n_^ 2 ) J ff(o)|o) (i) 

when ti > t 2 > 0. In cq. (§), IL/(i) is an inter- 
polating field of the final-state particle (denoted 
as "pion" in the following) with a fixed spatial 
momentum 



d 3 xe-^ s Il(x, t); 



(2) 



H (0) = H(x = 0, t = 0) is any local operator that 
couples to the two pions in the final state; T [. . .] 
represents the T-product of the fields and the vac- 
uum expectation value corresponds, in a numeri- 
cal simulation, to the average over the gauge field 
configurations. 
When t\ — * 00 

Gg(h,t 2 ) -> ^(o|n^i)|n)(n|n_ f (t 2 )ff(o)|o) 
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where 



E 9 =y/M* + &, V%i = (0|n(0)|g) 



(3) 
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and 

G 3 (t) = (q\n-&mm 



(5) 



for t > 0. 

Inserting a complete set of otif-states n, we can 
write 

(q\U-S)H(0)\0) =E( 2 ^ (3) (^) x 

n 

(q\IL(0)\n)(n\H(0)\0)e-( E "- E rt = 

^e- Eit out{q, -$|ff(0)|0) + ]T(27r) 3 X 
^ n 

S^(P n )(q\n(0)\n) c (n\H(0)\0)e-^- E ^(6) 



where the term proportional to ut(q, (0)\0) 
on the r.h.s. is the disconnected contribution. We 
now define the connected matrix element of the 
pion field as follows: 



(q\U(0)\n) c = 



2^[M(q,-q;n)}* 
{E n + 2Ef)(-E n + 2Ej-ie) 



(7) 



The LSZ reduction formula shows that, when 
the four momentum of the pion field 11(0) goes 
on the mass-shell (E n — > 2E^), M. reduces to 
the invariant scattering amplitude of the process 
n(q) + ir(—q) — > n, cf. eq. (15) of ref. S: 

(2n) i 5^(P n )5(E n - 2Ef)[iM(q, n)]* = 
in(q, -q\n)out -out (q, -q\n) ou t ■ (8) 



In general one could write 



(n\IL(0)\q) = ^Z n M% -q; n)T ( ^- 



(E + 2E ? - ie){~E + 2E ? - ie) 



(9) 



with the condition that the modulating factor T 
satisfies .F(l) = 1 for the on-shell pion (E n /2E^ = 
I). The factor in square brackets is (up to a fac- 
tor l/2Eg) the propagator of two non- interacting 
pions Q. 

Inserting the definition (Q) into eq. (||) , and us- 
ing the identity 



1 



E - 2E 3 - ie 
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we obtain 



,Egl 
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X out (q, -q\H(0)\0) + m (q, -q\H (0)\0) 



with 

ps) = (^)vj2< 



-(E n -2E s )t 



(27r) 3 ^(F„) 



M*(q,-q;n)(n\H(0)\0) 



(11) 



which corresponds to eq. (21) of ref . 0. Pq-(t) is 
a sum over off-shell amplitudes. In the limit t — > 
+oo, it is dominated by intermediate states with 
energy i£ n < 2_Eg-, which correspond to two-pion 
states with momenta k and —k, with |fc| < 
We therefore obtain, for t > 0, 



^(*) 



d 3 ki 



d 3 kn 



x (2nfd^(k 1 + k 2 )e-^ +E ^- 2E ^ 
M*{q,-q-,k u k 2 ){k u k 2 \H{Q)\(S) 
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(2tt)32^ (2^2^ 

x Ar(ft-$ifa,fe)(&,te|JT(o)|o> = 
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x (fc,-fc|iJ(0)|0) 



.M*(<f,-<f;/c,-fc) 
, • (12) 

\k\ = y?E 2 -4M%/2 



In the last equality, we have only considered 
s-wave two-pion scattering, and therefore inte- 
grated over the angles in the phase spacef]. 

3 In general, one should also include partial amplitudes 
with non-zero angular momentum, and retain the angu- 



Watson theorem states that, in absence of CP 
violation, 



out (k,-k\H(0)\0)=e iS ^A(s) 
and 

in (k,-k\H(0)\0) =e~ lS ^A(s) 



(13) 



(14) 



with A(s) real. The phase S is related to the two- 
pion scattering amplitude by 



(q,-$q,-q)i 



_ 2id(s 



(15) 



In the case of on-shell s-wave n — tt scattering, 
the invariant matrix element is given in terms of 
the phase by 



M(s) = 



i<5(a) 



sin<5(s) 



(16) 



We now make use of our "smoothness" hypothesis 
to write the off-shell matrix element as 

M(q, -q; k, -k) = {mn) ^ e^ smS(s) , (17) 



with s = 4.E|. Inserting the definitions 
fll?! ) into eqs. ([ll]) and (|l2|), we obtain 



and 
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G 3 (t) 



2M, 
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~ 4B| - E 2 



A(E 2 )smS(E 2 ) 



2E ? 



A(s q ) cos 6(s q ) + (4%) x 

dEe -{E-2E s )t ■ 

■A(s) sm S(s) 



2M, 



7T 4EI-E 2 



(18) 



with s q = AEl and s = E 2 . 



3. MESON DECAYS ON THE LATTICE 

To fix our ideas, let us consider D decays. In 
this case, the starting point is the four-point cor- 
relation function 

G(ti,t 2 ,t D ,q,-q,p D = 0) = 
= (U^U^Hw^dIMd)) (19) 



lar dependence. The generalization of the above formulae 
to the case of total momentum and angular momentum 
different from zero is straightforward. 



in the limit to — > — oo, t\ — > +00. Here Ti\y is the 
weak Hamiltonian. In this case eq. (|l8| ) becomes 

G(t 1 ,t 2 ,t D ,q,-q,p D = 0) = J^2^^ x 
e - E f(ti+ t 2)-M D |t D | |a w (s 8 ) cos + (4£ ? ) 



y^(s) sin<5(s) 



For numerical applications, it is convenient to 
consider the amputated correlation function given 
by the ratio 



R nw (t 2 ,q) 

where 

Sn (ti , Eq) - 



G(t 1 ,t 2 ,t D ,q,-q,p D = 0) 
S n (ti,E ? )S n (t2,E ? )S D {t D ,M D ) ' 



2E S 



(20) 



and similarly for the other meson propagators; 
R Hw (t 2 ,q) = S^A W \s q ) cos S(s q ) + (4£ f ) x 
r+0 ° dE e -(^-2%)t2 



2Mt tt 4£|-E2 



.4^(5) sin <5(s) 



(21) 



On a lattice of volume L 3 x T, eq. ( pl| ) becomes 
R Hw {t 2 ,q) = A(s q ) cos 5{s q ) + x 



E 



r {E,-2E s )t 2 



-A(s) sm6(s) 



, (22) 



where all the quantities are given in units of the 
lattice spacing, 

2tt 

Ei — \fs — 2Ej: with k = —(n x , n y ,n z ) ; (23) 

s q has been defined before; n XtViZ = 0, 1, . . . ,L — 1 
and J^g. denotes the sum over all the values 

of the energy corresponding to the momenta k 
allowed by the discretization of the space-time 
on a finite volume, excluding those correspond- 
ing to Ei = 2Eq-. Different combinations of mo- 
menta corresponding to the same energy should 



be included only once in the sum appearing in 
eq. (f22|). AEi = Ei + i — Ei is the difference 
of the nearest successive allowed values of Ei 
(Eq = 2M W , Eh = 2y 'Ml + (27r/La) 2 , E 2 = 
2^/M 2 + 2(2n/La) 2 , etc.). One can show that 
the expression in eq. ( p2| ) tends to the correspond- 
ing continuum one in eq. ( ^lf ) as L — > 00. How- 
ever, for presently available lattices, the allowed 
range of the momenta and of s q is limited, and 
therefore eq. (|22h might not be a good approxi- 
mation of eq. (|2l|). One might try to improve the 
approximation by using pn 7^ to increase the 
allowed values of the momenta, and by correct- 
ing for the difference between the discrete phase 
space and the continuum one. 

We now explain the strategy to extract the ma- 
trix element of Ttw from R Uw (t 2 ,q). To fix our 
ideas, we consider the case po = 0. The proce- 
dure is the following: 

1. Compute R nw (t 2 ,q = 0). The phase van- 
ishes at threshold, and therefore we have 



R nw (t 2 ,q = 0) = A w (4M£), 
and we can extract A w (AM 2 ). 



(24) 



2. Compute R Hw (t 2 ,q\) for the first allowed 
non-zero value of the momentum, \qi\ — 
(2ir) /L. The only value of k that might con- 
tribute in the sum in eq. (|22]), apart from 
exponentially suppressed terms, is k = 0, 
but at threshold the phase vanishes and 
therefore there is no contribution from the 
sum over k. Thus we have 

R Hw (t 2 ,qi) = A W (4E 2 ) cosS(AEl) , (25) 



with Ex = y/M* + (2ir) 2 /(La) 2 . 

3. Compute R Hw (t 2 ,q 2 ), with \q 2 \ = 
V2(2tt)/L. In this case, the term in the 
sum over k corresponding to \k\ = (2n)/L 
gives an exponentially increasing contribu- 
tion, and we get 

R Hw (t2,q 2 ) = A w (4El) cos 5(4El) + ^ x 

TT 

e 2(E 2 -E 1 )t 2 



(E 2 -E l )- 



A w (iE'f) sin S(4E(), 



Ei - Ef 



where E 2 = y/M% + 2(2ir) 2 / (La) 2 . 



4. CONCLUSIONS 



It is straightforward to derive the expression of 
R Uw for the next steps. In this way, we can ex- 
tract 



tan<5(s) 



As 



and \A(s)\ = VAj+A 2 (26) 



as a function of the centre-of-mass energy, where 
A c = A(s) cos 8(s) and A s = A(s) sin S(s) . 

An interesting case, which has been discussed 
in detail in ref. ||, is the one of FSFs domi- 
nated by the exchange of a resonance a in the 
s-channel^]. In this case, we can express the ma- 
trix elements of H and II in terms of the mass 
of the resonance M a , its width T = X/M a and 
its couplings g and V, to H and to the two-pion 
state respectively: 



out (q,~q\H (0)\0) 



(g|n(0)|fc,-fc) c 



Ml 



ix( Sq y 



\V(s)\< 



Ml 



iX(s) 



(E- 



2E r 



ie)(-E + 2E r 



(27) 



The phase shift and the amplitude are given in 
terms of the parameters of the resonance by the 
following relations: 



Ml - s 
and 

cosi5(s) 
sin<5(s) 



iX(s) 



A(s)e 



Ml 



iS(s) 



y/(Ml- S ) 2 +X( S ) 2 

X(s) 

y/(Ml- S ) 2 +X(s) 2 



(28) 

(29) 
(30) 



By studying the three-point H — II — II correla- 
tor Gq*(ti,t2), we can extract the parameters of 
the resonance and exploit this information when 
analyzing the four-point function. 

4 In ref. the result corresponding to eq. of the 

present work has been derived in the case of FSI's domi- 
nated by the exchange of a resonance. However, as we have 
shown above, the result ( |li| ) holds under a more general 
"smoothness" hypothesis. 



We have shown that, in spite of the MTNGT, it 
is in principle possible to extract the amplitude 
and phase of two-body non-leptonic decay ma- 
trix elements, under a reasonable "smoothness" 
hypothesis. We have sketched the strategy to 
extract the relevant information from Euclidean 
correlators computed numerically in lattice QCD. 
In the case of FSI's dominated by a resonance, it 
is possible to extract the parameters of the reso- 
nance from the study of a three-point correlation 
function. A feasibility study of the method we 
propose is currently under way on the APE ma- 
chine. 
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